Resummation for the Tevatron and LHC boson production at small x 



Stefan Berge 1 , Pavel Nadolsky 1 , Fredrick Olness 1 , and C.-P. Yuan 2 
Southern Methodist University, Department of Physics, Dallas, TX 75275-0175, U.S.A. 
2 Michigan State University, Dept. of Physics and Astronomy, East Lansing, MI 48824-1 1 16, U.S.A. 

Abstract 

Analysis of small-x semi-inclusive DIS hadroproduction suggests that multi- 
ple parton radiation leads to a broadening of transverse momentum (qx) distri- 
butions beyond that predicted by a straightforward calculation in the Collins- 
Soper-Sterman formalism. This effect can be modeled by a modification of the 
CN . resummed form factor in the small- x region. We discuss the impact of such 

C ! a modification on the production of electroweak bosons at hadron-hadron col- 

liders. We show that if substantial small-x broadening is observed in forward 
Z° boson production in the Tevatron Run-2, it would strongly affect predicted 
qx distributions for W, Z, and Higgs boson production at the Large Hadron 
Collider. 



o 



OS 



(N 



> 

In the production of electroweak bosons, precise knowledge of the transverse mass Mx and trans- 
verse momentum qx provides detailed information about the production process, including the mass of 
the boson and associated radiative corrections. At the Tevatron, qx distributions of Z° bosons offer in- 
sight into soft gluon radiation, and this information is then used for precision extraction of the W boson 
mass. At the LHC, good knowledge of the transverse distribution of Higgs bosons H° will be needed to 
efficiently separate Higgs boson candidates from the large QCD background. Accurate predictions for 
the small-c/T region are obtained via resummation of large logarithms \n n (qx/Q) arising from unsup- 
^ . pressed soft and collinear radiation in higher orders of perturbation theory. 

As we move from the 2 TeV Tevatron to the 14 TeV LHC, typical values of partonic momentum 
fractions x for producing W, Z°, and H° bosons become smaller, thus enhancing ln(l/x) terms in 
higher orders of a s . It is not entirely known how these terms (not included in a fixed-order cross section 
^ . or conventional qx resummation) will affect W, Z , and H° production at the LHC energies, in part 

because no Drell-Yan qx data is available yet in the relevant region of x of a few 10~ 3 or less. 

Studies Q |3 in the crossed channel of semi-inclusive deep-inelastic scattering (SIDIS) suggest 
that hadronic qx distributions at small x cannot be straightforwardly described within the Collins-Soper- 
Sterman (CSS) resummation framework 0, if the nonperturbative Sudakov function behaves like its 
large-x counterpart from the Drell-Yan process. A qx distribution in SIDIS at x < 1CP 2 is substantially 
broader than the conventional CSS prediction. The broadening effect can be modeled by including an 
extra x-dependent term in the Sudakov exponent. To describe the data, the extra term must grow quickly 
as x — > 0. It noticeably contributes to the resummed form factor at intermediate impact parameters (b ~ 
1/qx < 1 GeV -1 ), which hints at its origin from perturbative physics. A possible interpretation of this 
term is that it mimics higher-order contributions of the form a™ ln n (l/x), which are not included in the 
resummed cross section. Due to the two-scale nature of the qx resummation problem, the non-resummed 
ln(l/x) terms may affect the qx distribution even when they leave no discernible trace in inclusive DIS 
structure functions. The DIS structure functions depend on one hard scale (of order Q), while the CSS 
resummation formula (cf. Eq. Q) also includes contributions from large impact parameters b (small 
momentum scales). As b becomes large, the series a™(l/b) ln n (l/x) in the CSS formula may begin to 
diverge at a larger value of x than the series oOf (Q) ln n (1 / x) in the inclusive structure functions. For this 
reason, transition to /c^-unordered (BFKL-like [4|) physics may happen at larger x in qx distributions 
than in inclusive (one-scale) observables. 



Contributed to Workshop on Physics at TeV Colliders, Les Houches, France, 26 May - 6 June 2003. 



The qx broadening discussed above was observed in semi-inclusive DIS processes. In this study, 
we explore its possible implications for the (crossed) Drell-Yan process. We begin by examining the 
resummed transverse momentum distribution for the Drell-Yan process 0, following notations from 
Ref. 0: 

d^qj = If / Wf ^ Q ' XA ' XB) + Y{qT ' Q ' XA,XB) - (1) 

Here x a ,b = Qe ±y / ^/S, the integral is the Fourier transform of a resummed form factor W given in 
impact parameter (b) space, and Y is a regular (finite at qx — > 0) part of the fixed-order cross section. In 
the small-6 limit, the form factor W is given by a product of a perturbative Sudakov exponent e~ Sp and 
generalized parton distributions V(x,b): 



W(b,Q,x A ,x B ) =e- Sp ^V(x A ,b)V(x B ,b). (2) 



At moderately small x, where the representation © for W holds, we write these generalized parton 
distributions in the form 

^M)U<A- 2 ^(C®f)(x,b /b)e-P^ b \ (3) 

where C(x,bo/b) are coefficient functions, f{x,p) are conventional parton distributions, and bo = 
2 e -7B = 1.12... is a commonly appearing constant factor. 

The term e~ p ^ b2 in V(x,b) will provide an additional qx broadening, with an x dependence 
specified by p(x). For example, it may approximate x-dependent higher-order contributions that are not 
included in the finite-order expression for (C(g)/). We parametrize p(x) in the following functional form: 



p(x) = co L/-2+-2 > ( 4 ) 




such that p(x) ~ cq/x for x <C xq, and p(x) ~ for x ^> xo- This parameterization ensures that the 
formalism reduces to the usual CSS form for large x (x 3> a?o) an d introduces an additional source of qx 
broadening (growing as 1/x) at small x (x <C ^o). The parameter Co determines the magnitude of the 
broadening for a given x, while xo specifies the value of x below which the broadening effects become 
important. In principle, Co and xo may depend on the hard scale Q; in this first study, we neglect this 
dependence. Based on the observed dependence p{x) ~ 0.013/x at x < 10~ 2 in SIDIS energy flow 
data [2], we choose cq = 0.013 and xq = 0.005 as a representative choice for our plots. 

As x — > 0, the additional broadening term in Eq. © affects the form factor W both at perturbative 
(6 < 1 GeV -1 ) and nonperturbative (b> 1 GeV -1 ) impact parameters. In addition, the resummed cross 
section contains conventional non-perturbative contributions from power corrections, which become im- 
portant at large impact parameters (b > 1 GeV -1 ). We introduce these corrections by replacing the 



impact parameter b in functions Sp and (C f) with a variable 6* = b/yl + b 2 / (0.25 GeV 2 ) 
and including a nonperturbative Sudakov exponent exp {— S^pib, Q)}. The function Swp(b,Q) is 
parametrized by a 3-parameter Gaussian form from a recent global fit to low-energy Drell-Yan and Teva- 
tron Run-1 Z° data [5 1. Combining all the terms, we have: 

<la ; e-&*(C®f)(x A ,bo/b.)(C®f)(x B ,bo/h) 



dydq T S J (2tt) 2 

-Sp(b,,Q)-S NP {b,Q)-b 2 p(x A )-b 2 p(x B ) _|_ y ^ 



Figs.^and|2]show the comparison of the resummed cross section © with the additional broaden- 
ing term (p(x) ^ 0) to the resummed cross section without such a term (p(x) = 0). We consider cross 



pp -» Z°X (Vs= 1960 GeV) pp -> Z° X (VsT = 1 960 Ge V) 
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Fig. 1: qT distributions of Z° bosons in the Tevatron Run-2; (a) integrated over the full range of Z boson rapidities; (b) in- 
tegrated over the forward regions \y\ > 2. The solid curve is a standard CSS cross section, calculated using the 3-parameter 
Gaussian parametrization (5) of the nonperturbative Sudakov factor. The dashed curve includes additional terms responsible 
for the qT broadening in the small-a; region (cf. Eq. J5J). 



sections for the production of Z° and H° bosons, calculated according to the procedures in Refs. [6| and 
0, respectively. The numerical calculation was realized using the programs Legacy and ResBos OED, 
and with the CTEQ6M1 parton distribution functions [9|. The perturbative Sudakov factor was included 
up to 0{a 2 s ), and the functions (C f ) up to 0(a s ). The relevant perturbative coefficients can be found 
in Refs. EH. 

Fig- 02 a) shows the differential distribution da / dqx for Z boson production in the Tevatron Run-2, 
integrated over the rapidity y of the Z bosons. We observe that the distribution © with the additional 
small-x term (the dashed curve) essentially coincides with the standard CSS distribution (the solid curve). 
When y is integrated over the full range, both resummed cross sections are dominated by contributions 
from x ~ 0.05 3> Xq, where the additional broadening (given by the function p(x)) is negligible. For 
this reason, the Tevatron distributions that are inclusive in y (e.g., the Run-1 Z° boson data) will not 
be able to distinguish the small- x broadening effects from uncertainties in the nonperturbative Sudakov 
function Snp- 

In contrast, the small-x broadening does lead to observable differences in the distributions in 
the forward rapidity region, where one of the initial-state partons carries a smaller momentum fraction 
than in the central region. Fig. ^b) shows the cross section da/dqx for Z bosons satisfying \y\ > 2. 
The peak of the curve with p(x) ^ is lower and shifted toward higher qr. While this difference was 
not large enough to be observed in the Tevatron Run- 1 , it seems to be measurable in the Run-2 given the 
improved acceptance and higher luminosity of the upgraded Tevatron collider. The small-a; broadening 
is more pronounced in W boson production due to the smaller mass of the W boson. 

We now turn to the LHC, where the small-x broadening may be observed in the whole rapidity 
range due to the increased center-of-mass energy. Fig. |2a) displays the distribution da/dqx for Z° 
production with and without the small-x terms. Here, the difference is striking even if y is integrated out. 
Effects of a similar magnitude are present in W boson production, and they are further enhanced in the 
forward regions. 

The small-x broadening is less spectacular, but visible, in the production of light Standard Model 




Higgs bosons via the effective ggH vertex in the limit of a heavy top quark mass. Fig.|2b) displays the 
resummed cross sections for production of Higgs bosons with a mass Mh = 120 GeV for several choices 
of Snp and the broadening term. We first compare cross sections for p(x) = and p(x) ^ (thick 
lines), where the functions p(x) and SNp(b, Q) are taken to be the same as in Z° boson production. The 
difference between the two cross sections is not large, due to a harder qp spectrum in the Higgs boson 
case. The peaking in the gg-dominated H° distribution occurs at qp = 10 — 20 GeV, i.e., beyond the 
region where the function p(x)b 2 play its dominant role. This is different from the gg-dominated Z° 
distribution, where the peak is located at qp ~ 5 — 10 GeV and is strongly affected by exp {— p(x)b 2 }. 
Hence, for the same function p(x) as in the Z° boson case, the difference between the curves with and 
without p{x) is minimal. 

The harder qp spectrum in the Higgs boson case is induced by a larger leading-logarithm coef- 
ficient (Ca) in gg channels, as compared to the leading-logarithm coefficient Cp in qq channels. This 
suggests that the Q-dependent part (and possibly other terms) of the nonperturbative Sudakov function 
Snp in Higgs boson production is also multiplied by a larger color factor than in the Drell-Yan process. 
We estimate this effect by multiplying Sn p by the ratio of the leading color factors in Higgs and Z° bo- 
son production processes, Ca/Cf = 9/4 (the thin solid line). The resulting change turns out to be small 
because of the reduced sensitivity of the Higgs boson cross section to nonperturbative contributions. 

The ln(l/x) terms may be enhanced in the case of the Higgs bosons as well, due to the direct 
coupling of the Higgs bosons to gluon ladders. At present, we do not have a reliable estimate of the 
small-x broadening in gluon-dominated channels. However, this broadening would have to be quite 
large to affect qp of 10-20 GeV or more, i.e., in the region where selection cuts on the qp of the Higgs 
boson candidates will be imposed. For example, increasing the function p{x) by a factor of two as 
compared to the Z® boson case would lead to a distribution shown by the thin dashed line. While at 
qp > 20 GeV this effect is relatively small as compared to other theory uncertainties (e.g., higher-order 
corrections), it may affect precision calculations of qp distributions needed to separate the Higgs boson 
signal from the background in the 77 mode. 

Additional constraints on the small- x behavior of the resummed cross sections in the gg channel 
could be obtained from examination of photon pair production away from the Higgs signal region. As 



the mass of the photon pair decreases, 77 production in the gluon fusion channel via a quark box diagram 
becomes increasingly important. For instance, the subprocess gg — ► 77 contributes up to 40% of the 
total cross section at Q = 80 GeV flOl . By comparing qx distributions in pp — ► 77 and pp —> Z in the 
same region of Q, one may be able to separate the qq and gg components of the resummed cross section 
and learn about the x dependence in the gg channel. 

To summarize, we argue that a measurement of transverse momentum distributions of forward Z 
bosons at the Tevatron will provide important clues about the physics of QCD factorization and possibly 
discover broadening of qT distributions associated with the transition to small- x hadronic dynamics. 
Based upon the analysis of qx broadening effects observed in semi-inclusive DIS, we have estimated 
similar effects in the (crossed) processes of electroweak boson production at hadron-hadron colliders. 
While the estimated impact on the Higgs boson cross section da / dqx at high qx was found to be minimal, 
much larger effects may occur in W and Z boson production in the forward region at the Tevatron Run-2, 
and at the LHC throughout the full rapidity range. If present, the small- x broadening will have to be taken 
into consideration in precision studies of electroweak boson production. Additionally, its observation will 
provide insights about the transition to fc^-unordered (BFKL-like) dynamics in multi-scale distributions 
at hadron-hadron colliders. 
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